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Abstract The quantum vacuum may become an electromagnetic superconductor in 
the presence of a strong external magnetic field of the order of 10 16 Tesla. The mag- 
netic field of the required strength (and even stronger) is expected to be generated for 
a short time in ultraperipheral collisions of heavy ions at the Large Hadron Collider. 
The superconducting properties of the new phase appear as a result of a magnetic- 
field-assisted condensation of quark-antiquark pairs with quantum numbers of elec- 
trically charged p mesons. We discuss similarities and differences between the 
suggested superconducting state of the quantum vacuum, a conventional supercon- 
ductivity and the Schwinger pair creation. We argue qualitatively and quantitatively 
why the superconducting state should be a natural ground state of the vacuum at the 
sufficiently strong magnetic field. We demonstrate the existence of the supercon- 
ducting phase using both the Nambu-Jona-Lasinio model and an effective bosonic 
model based on the vector meson dominance (the p-meson electrodynamics). We 
discuss various properties of the new phase such as absence of the Meissner effect, 
anisotropy of superconductivity, spatial inhomogeneity of ground state, emergence 
of a neutral superfiuid component in the ground state and presence of new topolog- 
ical vortices in the quark-antiquark condensates. 



To appear in Lect. Notes Phys. "Strongly interacting matter in magnetic fields" (Springer), edited 
by D. Kharzeev, K. Landsteiner, A. Schmitt, H.-U. Yee 



1 Introduction 



Quantum Chromodynamics (QCD) exhibits many remarkable properties in the pres- 
ence of a very strong magnetic field. The external magnetic field affects dynamics 
of quarks because the quarks are electrically charged particles. As a result, the mag- 
netic field enhances the chiral symmetry breaking by increasing the value of the 
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(quark) chiral condensate [1 1. The change in the dynamics of quarks is also felt by 
the gluon sector of QCD because the quarks are coupled to the gluons. Therefore, the 
magnetic field may affect the whole strongly interacting sector and influence very 
intrinsic properties of QCD such as, for example, the confinement of color 12 O. 

In order to make a noticeable influence on the strongly interacting sector, the 
strength of the magnetic field should be of the order of a typical QCD mass scale, 
eB ~ nix, where m K ss 140MeV is the pion mass. The corresponding magnetic 
field strength, B~3x 10 14 T, is enormous from a human perspective (IT = 10 4 G). 
However, such strong magnetic field can be achieved in noncentral heavy-ion colli- 
sions at Relativistic Heavy-Ion Collider (RHIC) |4j. At higher energies of the Large 
Hadron Collider (LHC), noncentral heavy-ion collisions may generate even higher 
magnetic field of eB ~ 15m^ (B ~ 5 x 10 15 T) |4j. And in ultraperipheral colli- 
sions - when two nuclei pass near each other without a real collision - the mag- 
netic field strength may even reach eB ~ (60. . . 100)ra 2 or, in conventional units, 
B-(2...3)xl0 16 T, Ref. 0. 

Despite the magnetic field is generated in a heavy-ion collision for a very short 
time, it may have observable consequences. In noncentral collisions, the magnetic 
field is generated together with a hot expanding fireball of quark-gluon plasma. 
Topological QCD transitions may lead to a chiral imbalance of the plasma, and, 
in turn, the chirally-imbalanced matter may produce an observable electric current 
along the axis of the magnetic field [6| driven by the chiral magnetic effect J7). 

In a finite-density (quark) matter the magnetic catalysis [ 1 1 may be reversed [8| 
and the phase diagram may be modified substantially 1 8 , 9 ]. In the absence of matter 
(i.e., in the vacuum), the external magnetic field affects the finite-temperature phase 
structure of the theory by shifting the critical temperatures and altering the strength 
of the confinement-deconfinement and chiral transitions (HO. 

The vacuum may also spontaneously become an electromagnetic superconductor 
if the magnetic field strength exceeds the following critical value iKflfTTI : 

B c ~10 16 Tesla or eB c ~ 0.6GeV 2 . (1) 

The counterintuitive superconductivity of, basically, empty space, should always be 
accompanied by a superfluid component lfl2l[T3l . We discuss these effects below. 

In Section|2]we describe the mechanism and the basic features of the vacuum su- 
perconductivity in a very qualitative way. We compare in details the vacuum super- 
conductivity with an ordinary superconductivity. We also highlight certain similari- 
ties of this exotic vacuum phase with a magnetic-field-assisted "reentrant supercon- 
ductivity" in condensed matter and the electric-field-induced Schwinger electron- 
positron pair production in the vacuum of Quantum Electrodynamics. 

In Section [3] the emergence of the superconducting phase is explicitly demon- 
strated both in a bosonic p-meson electrodynamics [ 14 1 and in an extended Nambu- 
Jona-Lasinio model (131 . Various properties of the superconducting state are sum- 
marized in the last Section. 
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2 Conventional superconductivity, vacuum superconductivity 
and Schwinger pair creation: differences and similarities 

2.1 Conventional superconductivity via formation of Cooper pairs 

Before going into details of the magnetic-field-induced vacuum superconductivity 
let us discuss certain very basic qualitative features of the conventional supercon- 
ductivity. Why certain compounds are superconductors? 

In a simplified pictur^] electrons in a metal can be considered as (almost free) 
negatively charged particles which move through a periodically structured back- 
ground of a lattice made of positively charged ions. The individual electrons scatter 
inelastically off the ions leading to a dissipation of an electric current and, conse- 
quently, to emergence of a nonvanishing electrical resistance in the metal. 

As an electron moves through the ionic lattice it attracts neighboring ions via a 
Coulomb interaction. The attraction leads to a local deformation of the ionic lattice, 
and, simultaneously, to an excess of the positive electric charge in a vicinity of 
the electron. The excess of the positive charge, in turn, attracts another electron 
nearby, so that in a background of the positively charged ion lattice the like-charged 
electrons may experience a mutual attractive force, Fig.[TJleft). The deformation of 
the ionic lattice can be viewed as a superposition of collective excitations of the ion 
lattice (phonons), so that the process of the electron-electron interaction via lattice 
deformations can be described by a phonon exchange. 

The attractive force between the electrons may, in principle, lead to formation 
of electron-electron bound states. However, this attraction is extremely weak and 




Fig. 1 (left) Formation of the Cooper pair (the yellowish oval) of electrons (the small green circles) 
in an ionic lattice (the large red circles) due to phonon interaction, (right) Two interacting electrons 
(the small green circles) and the Fermi sphere (the large blue circle) in the momentum space. The 
electrons in the Cooper pair have mutually opposite spins (the green arrows) and momenta. 



1 A comprehensive introduction to superconductivity can be found in books I16II17|[T81 . 
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therefore thermal fluctuations at room temperature easily destroy the two-electron 
bound states. On the other hand, at low temperature the attractive interaction be- 
tween the electrons prevails the thermal disorder and, consequently, the bound states 
may indeed be formed. These bound state are called the Cooper pairs. The electrons 
in the Cooper pairs have mutually opposite spins and opposite momenta thus mak- 
ing the Cooper pair a (composite) spin-zero bosonic state. Bosons have a tendency 
to condense at low temperature so that the Cooper pair condensate may emerge. 

In the condensed state all Cooper pairs behave as one collective entity. The 
Cooper-pair condensate can move frictionlessly through the ion lattice similarly 
to a motion of a superfluid. The motion without dissipation is guaranteed due to an 
energy gap, which separates the energy of the condensed ground state and a next, 
excited state. Since the intermediate states are absent, at certain conditions (low 
enough temperature, weak enough electric current, etc) the dissipative scattering of 
the Cooper pairs off the ions becomes kinematically impossible so that the motion of 
the Cooper pair condensate proceeds without dissipation. Since the Cooper pairs are 
electrically charged states, the condensation of the Cooper pairs turns the material 
into a superconducting state. 

The formation of the Cooper pairs is facilitated by the fact that at low tempera- 
ture the dynamics of the electrons becomes effectively one dimensional while in one 
spatial dimension even a very weak attraction between two particles should always 
lead to formation of a bound state (in the condensed matter context this property 
is known as the Cooper theorem). The effective dimensional reduction of the elec- 
tron dynamics from three spatial dimensions to one spatial dimension is possible 
because at low temperature the interaction between the electrons occurs if and only 
if the electrons are sufficiently close to the Fermi surface. One component of the 
momentum at the Fermi surface counts the degeneracy of the electron states while 
the other component is a dynamical degree of freedom. The Cooper pair is formed 
by two electrons with mutually opposite momenta and mutually opposite spins, Fig- 
ure [TJright). 

Summarizing, in order to exhibit the conventional superconductivity a system 
should satisfy the following basic requirements: 

A) electric charge carriers should be present in the system (otherwise the system 
cannot support the electric current); 

B) dynamics of the electric charge carriers should effectively be one-dimensional 
(otherwise the Cooper pairs cannot be formed); 

C) the like-charged carriers should experience mutual (pairwise) attraction (other- 
wise the Cooper pairs cannot be formed). 

Surprisingly, the same requirements are satisfied by the quantum fluctuations 
of the vacuum in a background of a sufficiently strong magnetic field. In the next 
section we compare basic features of conventional and "vacuum" superconductors. 
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2.2.1 Condition A: presence of electric charges 

In order for the vacuum to behave as an electromagnetic superconductor, one needs, 
at least, the presence of electrically charged particles (condition A on pageQ. From 
a first sight, it is impossible to satisfy this requirement because under the usual 
conditions the vacuum is characterized by the absence of the free electric charges. 
Nevertheless, the quantum vacuum may be considered as an excellent "reservoir" 
of various particles including the electrically charged ones. Moreover, under certain 
external conditions the virtual particles may become real. 

This "virtual-to-real" scenario does not sound unlikely. For example, there are 
at least two well-known cases of external conditions when a vacuum becomes an 
electrically conducting media: the vacuum may conduct electricity if it is either 
subjected under a strong electric field or if it is sufficiently hot. 

The first example of the "virtual-to-real" transition is the Schwinger effect in 
Quantum Electrodynamics (QED): a sufficiently strong external electric field gen- 
erates electron-positron pairs out of the vacuum |[T9l . The created positrons and 
electrons move in opposite directions thus creating an electric currenrl The critical 
strength of the electric field required for this process is E c = rn^/e « 10 18 V/m. 

The second "virtual-to-real" example is a simple thermal ionization of electron- 
positron pairs: the vacuum turns into an electron-positron plasma at temperatures 
T ~ 0.1 T QED where r QED « 2m e « 1 MeV « 10 10 K is a typical QED temperature. 

Table 1 Conventional superconductivity vs. vacuum superconductivity: very general features 
(from Ref. l20l ). 



Property 


Conventional superconductivity 


Vacuum superconductivity 


Environment 


a material (metal, alloy etc) 


vacuum (empty space) 


Reservoir of carriers 


real particles 


virtual particles 


Normal state 


a conductor 


an insulator 


Basic carriers of 
electric charge 


electrons (e) 


light quarks (u, d) and 
light antiquaries (u, d) 


Electric charges 
of basic carries 


q e = -e (e=|e|) 


q u = +2e/3, qa = -e/3 
q r , = -2e/3, qj = +e/3 



Thus, the quantum vacuum may be turned into a conductor if it is subjected to 
sufficiently strong electric field (E ~ 10 18 V/m) or to sufficiently high temperature 
(T ~ 10 9 K). Below we show that sufficiently strong magnetic field (B ~ 10 16 T) 

2 We briefly discuss an analogy between the magnetic-field-induced vacuum superconductivity and 
the Schwinger effect in Section [2.2.5| page 1 1 1 1 
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may turn the vacuum into a superconducting state. The magnetic-field-induced 
vacuum superconductivity works at the QCD scale: the key role here is played by 
virtual quarks and antiquarks which have fractional electric charges. As we discuss 
below, the strong magnetic field catalyses the formation of the electrically charged 
condensates made of quarks and antiquarks. Very general features of a conventional 
superconductor and the magnetic-field-induced vacuum superconductivity are sum- 
marized in Table Q] 



2.2.2 Conditions B and C: formation of superconducting carriers 

In order for the superconducting earners to be formed, the fermion dynamics should 
be reduced from three spatial dimensions to one spatial dimension (condition B of 
Section 2.2.1 page|4|i. In conventional superconductivity the dimensional reduction 



proceeds via formation of the Fermi surface at sufficiently low temperatures. This 
mechanism cannot work in our case because the Fermi surface, obviously, does not 
exist in the vacuum due the very absence of matter. However, the dimensional reduc- 
tion may be achieved with the help of a magnetic field background since electrically 
charged particles with low-energy can move only along the axis of the magnetic 
field. This effect leads to the required dimensional reduction of the charge's dynam- 
ics from three to one spatial dimensions. 

The described dimensional reduction effect in the background of the external 
magnetic field works for all electrically charged elementary particles, including 
electrons, positrons, quarks, antiquarks etc. However, the superconducting bound 
state may only be formed from a particular combination of these particles which 
should satisfy the following conditions: 

(i) the superconducting bound state should be a boson; 

(ii) the bound state should be electrically charged; 

(iii) the interaction between the constituents of the bound state should be attractive. 

Condition (i) implies that the superconducting bound state should contain even 
number of constituents because the known carriers of the electric charge are fermions 
(quarks, electron and positron, etc). Below we consider simplest, two-fermion states. 

Condition (ii) implies that the bound state cannot be composed of a particle and 
its antiparticle. In combination with condition (iii) it means that the vacuum su- 
perconductivity cannot - unlike the Schwinger's pair creation - emerge in the pure 
QED vacuum sector which describes electrons, positrons and photons. Indeed, the 
electron-electron interaction is mediated by a repulsive photon exchange so that 
condition (iii) is not satisfied. On the other hand, the interaction between electron 
and positron is attractive, but the electron-positron bound state is electrically neutral 
so that in this case condition (ii) is not satisfied. Thus, the superconductivity cannot 
emerge in the pure QED. 

Therefore, the candidates for the superconducting charged bound states should 
be outside of the purely electrodynamics sector. Below we concentrate on the next 
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(in terms of energy scale), strongly interacting sector which describes the dynamics 
of quarks and gluons. 

The QCD sector of the vacuum contains the gluon particle which is a carrier 
of the strong force. From our perspective the gluon is an analogue of the phonon 
of conversional superconductivity because it is the gluon may which may provide 
an attractive interaction between quarks and antiquarks regardless of their electric 
charges. In particular, the gluon may bind a quark and an anti-quark into an elec- 
trically charged meson. The attractive nature of the gluon interaction allows us to 
satisfy condition C of superconductivity on page [4] 

Thus, the suggested mechanism of the vacuum superconductivity may indeed 
work at the interface of the QED and QCD sectors. The simplest example of the 
superconducting carrier may be given by a bound state of a u quark with the electric 
charge q u — +2e / 3 and a d antiquark with the electric charge q j = — q c ( = +e/3 . The 
attractive nature of the gluon-mediated interaction between the quark and antiquark 
of different flavors is only possible if these constituents reside in a triplet state, so 
that the ud bound state should be a spin-1 state (the p meson). 

Therefore, the vacuum analogue of the Cooper pair are the charged p* meson 
states. And in next sections we show that the p-meson condensates do indeed appear 
in the vacuum in the presence of the strong magnetic field, and we argue that the 
emergent state is indeed an electromagnetic superconductor. 



Table 2 Conventional superconductivity vs. vacuum superconductivity: superconducting carriers 



Property of carrier 


Conventional superconductivity 


Vacuum superconductivity 


Type 


Cooper pair 


p-meson excitations, p + and p~ 


Composition 


electron-electron state (ee) 


quark-antiquark states 
(p + = ud and p~ = du) 


Electric charge 


-2e 


+e and —e, respectively 


Spin 


typically spin-zero state (scalar) 


one-one state (vector) 




1) reduction of dynamics of basic electric charges 


The carriers are 


from three spatial dimensions to one dimension, 3d — > Id 


formed due to 


2) attraction force between 


2) attraction force between 




two electrons 


a quark and an antiquark 


1) a reason for the 
reduction 3>d — > Id 


at very low temperatures 
electrons interact with each 


in strong magnetic field the 
motion of electrically charged 




other near the Fermi surface 


particles is one dimensional 


2) attraction is due to 


phonons (lattice vibrations) 


gluons (strong force, QCD) 


Isotropy of 

superconducting 

properties 


Yes: superconducting 
in all spatial directions 


No: superconducting along 
the axis of the magnetic field, 
insulator in other directions 
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Notice that the formation of the bound state is facilitated by the dimensional 
reduction of the quark's dynamics in the background of the magnetic field (condi- 
tion B). The dimensional reduction implies automatically a strong anisotropy of the 
suggested superconductivity since the electric charges (the quarks u and d and their 
antiquarks) may move only along the axis of the magnetic field. As a result, the 
superconducting charge carriers (the p mesons in our case) may also flow along the 
axis of the magnetic field only. Thus, the vacuum exhibit a superconducting property 
in the longitudinal direction (along the magnetic field) while in the two transverse 
directions the superconductivity of the vacuum should be absent. 

Note that due to the anisotropic superconducting properties the vacuum in the 
strong magnetic field acquires a very unusual optical property: the vacuum becomes 
as (hyperbolic) metamaterial which behaves as diffractionless "perfect lenses" ETIl . 

In Table [2] we compare of certain basic features of the superconducting carriers 
in a conventional (low-temperature) superconductivity and in the vacuum (high- 
magnetic-field) superconductivity. 

2.2.3 Counterintuitive coexistence of magnetic field and superconductivity 
due to strong anisotropy of magnetic-field-induced superconductivity 

So far we have ignored a well-known property of all known superconductors: 

• Magnetic field (regardless of its strength) and conventional superconductivity^] 
(regardless of its mechanism) cannot coexist with each other! 

Thus, we can ask ourselves: why do we believe that the superconducting phase of 
the vacuum can exist in (and, moreover, be induced by) the strong magnetic field? 
In fact, this single question contains two puzzles (Table[3]l: 

• Why the Meissner effect is absent in the superconducting phase of vacuum? 

• Why strong magnetic field does not destroy the superconductivity of vacuum? 

A short "technical" answer to these questions is that in the background of the 
magnetic field the superconducting state of the vacuum has lower energy compared 
to the energy of the normal (insulator) state (Section |XT| page[T2]i. A physical argu- 
ment is that the strong magnetic field may coexist with the vacuum superconductiv- 
ity because the latter is highly anisotropic. Let us consider this point in detail. 

Qualitative arguments against the Meissner effect in the vacuum superconductor 
are as follows. The Meissner effect is a screening of weak external magnetic field by 
a superconducting state so that a magnetic field cannot penetrate deeply into a super- 
conductor. Qualitatively, the Meissner effect is caused by superconducting currents 
which are induced by the external magnetic field in the bulk of a superconductor. 
The circulation of these currents in the transversal (with respect to the magnetic 
field axis) plane generates a backreacting magnetic field, which screens the external 
magnetic field in the bulk of the superconducting material. The backreacting cur- 
rents are geometrically large, so that the corresponding magnetic length (i.e., the 



3 Except for the unconventional reentrant superconductivity (Section 2.2.6 page 
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Table 3 A comparison of the effects of magnetic field and thermal effects on conventional super- 
conductivities and electromagnetic superconductivity of vacuum (from Ref. fM ). 



Property 


Conventional superconductivity 


Vacuum superconductivity 


Magnetic field 


destroys superconductivity 


induces superconductivity 


The Meissner effect 


present 


absent 


Thermal fluctuations 


destroy superconductivity 


destroy superconductivity 



radius of the lowest Landau level), 1 / yjeBf, is much larger than the correlation 
length ^ of the superconductor. Since the vacuum superconductivity is realized only 
along the axis of the magnetic field, the large transversal currents are absent and the 
Meissner effect cannot be realized. 

If the axis of the external magnetic field is oriented along the normal to a bound- 
ary of an ordinary superconductor, then the backreacting magnetic field squeezes 
the external magnetic field into thin Abrikosov vortice^Jwhich form a sparse vortex 



lattice in a background of weak magnetic field. As we show below (Section 3.4.4 



page 25 1, the superconducting ground state of the vacuum is a dense lattice of 



Abrikosov-type vortices for which the magnetic length is of the order of (or even 
smaller than) the correlating length. This is a quantum regime of the Abrikosov 
lattice, in which the geometrically short transverse currents in the cannot screen ge- 
ometrically large external magnetic field. In this case the physical situation is similar 
to a "reentrant superconductivity"of extreme type-II superconductor in a high mag- 
netic field 123 (Section |ZZ6} pagefTl). 

If the axis of the weak external magnetic field is oriented tangentially to a bound- 
ary of an ordinary superconductor then the external field is usually expelled from 
the superconductor's bulk without formation of the Abrikosov vortices. In our case, 
the superconducting ground state is created by a strong magnetic field, and there- 
fore the very imposition of a weak tangential magnetic field is logically impossible 
from very simple geometrical reasons: a sole result of the superposition of the weak 
"testing" magnetic field onto the strong "creating" magnetic field is a slight turn of 
the stronger field. In other words, the weak testing field should slightly reorientate 
the anisotropy axis of the superconductor without destroying it. 

The ordinary superconductivity is destroyed by sufficiently strong magnetic field. 
Qualitatively, one can understand this effect as follows: in strong enough external 
magnetic field the (positive) excess in energy of the induced transverse supercon- 
ducting currents prevails the (negative) condensation energy of the superconducting 
carriers. As a result, at certain critical field the conventional superconductivity be- 
comes energetically unfavorable and the material turns from the superconducting 
state back to the normal (nonsuperconducting) state. On the contrary, in the vac- 
uum superconductivity the large superconducting currents are absent due to strong 



4 Since the magnetic flux coming through the superconductor's boundary is a conserved quantity, 
the superconductor expels it from the superconductor's bulk into thin vortexlike structures. 
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anisotropy of the superconducting currents, so that the mentioned argument should 
not work. Moreover, as we discuss below, the energy of the short transverse currents 
is diminished as the magnetic field becomes stronger. 

Thus, the electromagnetic superconductivity of the vacuum coexists with high 
magnetic field due to the anisotropy of the magnetic-field-induced superconducting 
properties. It is the anisotropy which makes the vacuum superconductivity to be 
different from the conventional one. 



2.2.4 Magnetic-field-induced vacuum superconductivity: temperature effects 

The common key element of the ordinary and vacuum superconductivities is the di- 
mensional reduction of the dynamics of the charge carriers (condition B on page[4]). 
Thermal fluctuations should destroy this property regardless of the mechanism of 
the dimensional reduction. In ordinary superconductivity, if the energy of the ther- 
mal fluctuations becomes of the order of the Fermi energy then the Fermi surface 
broadens and the dimensional reduction no more works. 

The thermal fluctuations should destroy the vacuum superconductivity because of 
the same reason as in the ordinary superconductivity: the loss of the dimensional re- 
duction (Table|3]l. Indeed, the one-dimensional motion of electric charges in strong 
magnetic field is realized due to the fact that that the electric charges occupy the 
lowest Landau level which is localized in the transverse plane. The one-dimensional 
motion can only be spoiled by transitions of the particles to higher, less localized 
Landau levels. Generally, for a typical gap between the Landau energy levels is ex- 
pected to be of the QCD scale, 8E ~ Aqcd ~ lOOMeV, so that thermal fluctuations 
of a typical QCD scale, T ~ Aqcd should destroy the dimensional reduction. 

Therefore, we conclude that the superconductivity should be lost at certain crit- 
ical temperature, T c = T C (B). At the critical magnetic field, B = B c , the critical 
temperature is zero, T C (B C ) = 0. The corresponding phase diagram is schematically 
shown in Fig. [2] the superconducting and hadronic phases are separated by a phase 
transition of (presumably) second order ifTTl . 




Hadronic phase 



* 



* 




electromagnetic 
superconductivity 



transition. 



on d order) 

Phase of 







B c 



Magnetic field, B 



Fig. 2 Schematic plot of the suggested QCD phase diagram in the presence of magnetic field in 
the low temperature region 1111 . 
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2.2.5 Electric-field-induced pair production (the Schwinger effect) and 
magnetic-field-induced superconductivity: a comparison 

The Schwinger effect is a generation of the electron-positron pairs from the vacuum 
in a background of a strong enough electric field [ 19 1. The created particles form a 
momentary electric current which tend to screen the external electric field which has 
created them. The electron-positron pair production is a process which is described 
entirely by the QED sector of the vacuum. 

The vacuum superconductivity is associated with the emergence of the electri- 
cally charged quark-antiquark condensates out of vacuum provided the vacuum is 
subjected to the strong enough magnetic field ifTOlfTTI . Contrary to the Schwinger 
effect, these electrically charged condensates do not screen the external magnetic 
field which has created them. 

Following Ref. [20 1, in Table [4] we compare the very basic features of the 
Schwinger effect and the vacuum superconductivity. 



Table 4 Basic features of the Schwinger effect and the electromagnetic vacuum superconductivity. 



Property 


Schwinger effect 


Vacuum superconductivity 


Environment 


vacuum 


vacuum 


Background of 


strong electric field, E 


strong magnetic field, B 


Interactions involved 


electromagnetic (QED) 
only 


electromagnetic (QED) 
and strong (QCD) 


Typical energy scales 


megaelectronvolts (10 6 eV) 


gigaelectronvolts (10 9 eV) 


Critical value 


E c = ml/e^ 10 18 V/m 


B c = m 2 p je^ 10 16 T 




(m e — 0.5 1 1 MeV is electron mass) 


(nip — 0.775 GeV is p-meson mass) 




virtual electron-positron (e~e + ) 


virtual quark-antiquark pairs 


Nature of the effect 


pop up from the vacuum 


(km and dd) pop up and form 




and become real e~e + pairs 


real ud and dii condensates 


Backreaction 


created e~e + pairs tend 


created ud and dii condensates 


to screen the external field 


do not screen the external field 


Stability 


a process (unstable) 


a ground state (stable) 


Transport property 


an electromagnetic conductor: 
electric current is generated 


a steady superconducting state 



2.2.6 Electromagnetic superconductivity of vacuum and "reentrant 
superconductivity" in strong magnetic field 

The magnetic-field-induced electromagnetic superconductivity of vacuum may have 
counterparts in certain condensed matter systems. It was suggested in Ref. Il23ll that 
in a very strong magnetic field the Abrikosov flux lattice of a type-II superconductor 
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may enter a quantum limit of the "low Landau level dominance", characterized by 
a spin-triplet pairing, absence of the Meissner effect, and a superconducting flow 
along the magnetic field axis. The mentioned quantum limit is reached when the 
magnetic length 1 / becomes of the order of the correlation length t, . 

In condensed matter, the magnetic-field-induced anisotropic superconductivity 
is sometimes called the "reentrant" superconductivity because the system should 
normally "exit" a superconducting state as an increasing external magnetic field 
suppresses superconductivity via diamagnetic and Pauli pair breaking effects. Al- 
though it is unclear whether this particular mechanism of the reentrant supercon- 
ductivity works in real superconductors, the magnetic-field-induced restoration of 
superconductivity was experimentally observed in certain materials like an uranium 
superconductor URhGe \22\. 

Our proposal ifTOl [TT1 of the vacuum superconductivity has basically the same 
features as the reentrant superconductivity 1231 : the electrically charged condensates 
correspond to a spin-one quark-antiquark states (p mesons), the vacuum supercon- 
ductor exhibits no Meissner effect while the vacuum superconductivity is highly 
anisotropic. 



3 Ground state of vacuum superconductor 

3.1 Energetic favorability of the superconducting state 

As we have argued in the previous section, a quark-antiquark pair of different flavors 
may condense in sufficiently strong magnetic field. How strong should the relevant 
magnetic field be? Following Ref. 11011 . let us make a simple estimation of the criti- 
cal magnetic field B c using very general arguments. 

An electromagnetic superconductivity emerges when an electrically charged field 
starts to condense. Assume that we have a free relativistic particle with mass m, 
electric charge e and spin s. In a background of a constant uniform magnetic field B 
the relativistic energy levels e of this particle are given by the following formula: 

«?U (P*) = Pz + ( 2n - + l)\eB\+ m 2 , (2) 

where n ^ is the nonnegative integer, s z = — 5, . . . ,s is the projection of the spin s 
on the field's axis, p z is the particle's momentum along the field's axis and g is the 
gyromagnetic ratio (or, "g-factor") of the particl^] 

Let us consider quark-antiquark bound states made of lightest, u and d quarks. 
The corresponding simplest spin-zero and spin-one bound states are called % and p 
mesons, respectively. 

5 Here we ignore the internal structure (formfactors) of the mesons treating them as pointlike bound 
states. We proceed similarly to the case of the ordinary superconductivity, where the large-sized 
Cooper pairs can also be treated as local objects in certain approaches, Section |3~2| page[l4] More 
general and formfactor-independent treatment is presented in Section|3.5| page|30| 
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The ground state energy (or, "mass") of the % mesons in the background of the 
magnetic field corresponds to the quantum numbers p z = and n z = (notice that 
s z = since the % meson is a spineless particle): 

ml ± {B)=ml + \eB\, (3) 

where m K = 139.6 MeV is the mass of the % meson in the absence of the magnetic 
field. The mass of the neutral n° meson is insensitive to the external magnetic field in 
our approximation, m K o (B) = m K (here we ignore a small splitting between masses 
of charged and neutral % mesons at B — 0). 

Analogously, the ground state energy ("mass") of the charged p ± mesons corre- 
spond to the quantum numbers p z = 0, n z = and s z = 1 : 

m 2 p± (B)=m 2 p -\eB\, (4) 

where m p = 775.5 MeV is the mass of the charged p meson in the absence of the 
magnetic field. The mass of the neutral p° meson is a B-independent quantity in this 
approximation, m p o (B) = m p (a small difference in masses of charged and neutral 
p mesons at B = is ignored again). 

It is important to mention that in Eq. |4]) the gyromagnetic ratio of the p meson 
is taken to be g — 2. This anomalously large value was independently obtained in 
the framework of the QCD sum rules [24] and in the Dyson-Schwinger approach 
to QCD 11251 . It was also conformed by the first-principle numerical simulations of 
lattice QCD ll26ll providing us with a value g ~ 2. The condensation of the charged p 
mesons in the vacuum of QCD is very similar to the Nielsen-Olesen instability of the 
pure gluonic vacuum in Yang-Mills theory [27 1 and to the magnetic-field-induced 
Ambj0rn-Olesen condensation of the W -bosons in the vacuum of standard elec- 
troweak model [28 1. Both the p mesons in QCD, the gluons in Yang-Mills theory, 
and the W bosons in the electroweak model have the anomalously large g-factor, 
g ~2. Notice, that the phase diagrams of QCD at finite density (at finite chemi- 
cal and/or isospin potential) contain certain phases characterized by the presence of 
exotic vector condensates ||9] [29] [30] [3T) . Some of these phases exhibit supercon- 
ducting/superfluid behavior [30, 31 1. 

In the absence of the magnetic field background the p meson is a very unstable 
particle. One can notice, however, that the ground-state mass Q of the charged n 
meson is an increasing function of the magnetic field strength B while the mass of 
the charged p mesons Q is a decreasing function of B. As all known modes of 
the p ± -meson decays proceed via emission of the ft* mesons, p* — > n ± X l32ll . 
it is clear that at certain strength of magnetic field the fast hadronic decays of the 
p mesons become forbidden due to simple kinematical arguments (the mass of the 
would-be decay products exceed the mass of the p meson itself). Thus, in a back- 
ground of sufficiently strong magnetic field the charged p meson should be stable 
against all known hadronic decay modes |[T0l . 

The presence of the superconducting ground state at high magnetic field can be 
seen as follows. The square of mass of the charged p meson should decrease as 
the magnetic field B increases, Eq. Q. When the magnetic field reaches the critical 
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value dl), the ground state energy of the p mesons becomes zero. As the magnetic 
field becomes even stronger, the ground state energy becomes a purely imaginary 
quantity indicating the presence of a tachyonic instability of vacuum. In other words, 
the trivial ground state, (p) =0, is no more stable at B > B c , and the vacuum should 
slide towards a new state with a nonzero p-meson condensate, (p) ^ 0. Since the 
condensation of the electrically charged field indicates the presence of electromag- 
netic superconductivity, the vacuum should become a superconductor at B > B c . 



3.2 Approaches: GinzJburg- Landau vs Bardeen-Cooper-Schrieffer 

The condensation of the p mesons in QCD in the background of the strong magnetic 
field may be treated in the same way as the condensation of the Cooper pairs in the 
conventional superconductivity. 

The conventional superconductivity may be described in the framework of both 
microscopic fermionic models and macroscopic bosonic theories. The fermionic 
models of the Bardeen-Cooper-Schrieffer (BCS) type describe basic carriers of elec- 
tric charge (electrons and/or holes). The BCS models are, generally, nonrenormal- 
izable because their Lagrangians include four-fermion term(s). The bosonic models 
of the Ginzburg-Landau (GL) type are usually based on renormalizable effective 
Lagrangians which describe superconducting excitations (the Cooper pairs) ifTTl . 

Despite of the fact that fermionic and bosonic approaches are formulated in a 
very different way, they both can describe the superconductivity at a good quan- 
titative level. Moreover, the bosonic and fermionic approaches are mathematically 
equivalent near the superconducting phase transition f33l . In Table [5] we outline a 
correspondence between the traditional (BCS and GL) models of conventional su- 
perconductivity and their vacuum counterparts which will be used to describe the 
magnetic-field-induced vacuum superconductor later. 



Table 5 Simplest models which are used to study physical properties of conventional and vacuum 
superconductivities . 



Basic field describes ... 


conventional 


vacuum 


bosonic, condensed 
superconducting carriers 


Ginzburg-Landau 
model t34l 


p-meson electrodynamics |14| based 
on vector dominance model 1351 


fermionic constituents of 
superconducting carriers 


Bardeen-Cooper- 
Schrieffer model 1371 


Nambu-Jona-Lasinio model 1361 
extended with vector interactions 1151 



It is important to notice that in the effective GL approach the Cooper pairs are 
treated as pointlike particles. However, physically the Cooper pairs are rather non- 
local objects because their size is much larger than the average distance between 
electrons in metals. Nevertheless, the GL model describes superconductivity very 
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well especially near the second-order phase transition where the symmetries of the 
system dominate its dynamics according to the universality argument. 

Following our experience in the conventional superconductivity, we describe be- 
low a GL like approach to the vacuum superconductivity using a model of a p- 
meson sector of vacuum [ 14 1. Then, we briefly outline a Bardeen-Cooper-Schrieffer 
approach to the vacuum superconductor using a well-known extension [15] of the 
Nambu-Jona-Lasinio model 1361 . We would like also to mention that signatures of 
the vacuum superconductor were also found in holographic effective theories 1 39 1 
and in numerical ("lattice") approaches to QCD BOl . 



3.3 Example: Ginzburg-Landau model 

Before going into the details of the p condensation in QCD, it is very useful to out- 
line a few basic properties of the conventional superconductivity in the GL model. 

3.3.1 The relativistic version of the Ginzburg-Landau Lagrangian 

Conventional superconductors can be described by the following GL Lagrangian: 

^gl = -^V v + (S)^)*2^-A(|0| 2 -t7 2 ) 2 , (5) 

where = — ieA^ is the covariant derivative and <P is the complex scalar field 
carrying the electric charg^Je. 

The superconducting ground state of the GL model, rj 2 > 0, is characterized 
by the presence of the homogeneous condensate (4>), with |(<£)| = rj, which is 
responsible for the superconductivity. The condensate breaks the electromagnetic 
symmetry group of Lagrangian <t> —> e' em <t> and —> + (0. 

In the superconducting phase the mass of the scalar excitation, 8<P = <P — (<P), 
and the mass of the photon field A^ are, respectively, as follows: 

m<p = V4Ar/ , m A = \Jler\ . (6) 

The classical equations of motion of the GL Lagrangian Q are: 

S) m £>' 1 + 2A(|0| 2 -tj 2 )0 = 0, (7) 

d v F v »+j£ L = 0, (8) 

where the electric current is 

J^ L = -ie[<P*D ll <P-(Ti^<P)*<P] . (9) 



6 Without loss of generality we consider the singly-charged bosons <P instead of the doubly charged 
Cooper pairs and we use a relativistic description of superconductivity. 
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Thermal fluctuations make the condensate (<f>) smaller, eventually destroying 
the superconductivity at certain critical temperature, T = T c , In order to describe 
this effect in the GL approach, one usually assumes that the quadratic coefficient of 
the potential term Vgl = •^(l^'l 2 — I ? 2 ) 2 m the GL Lagrangian ^ exhibits a linear 
temperature dependence: 

y GL (^) = ao(r-r c )|^| 2 +A|0| 4 + const. (10) 

The superconducting condensate is present at T < T c , while (<£>) = at T ^ T c . 



3.3.2 Magnetic field destroys conventional superconductivity 

In a background of a sufficiently strong magnetic field, B > B^ L , the superconduct- 
ing condensate disappears. The corresponding critical value of the magnetic field, 

B? L = ^-rj 2 , (11) 
le e 

can be obtained as follows. Consider a near-critical case when the uniform time- 



independent magnetic field B = F\i is slightly weaker than the critical value ([TT 
< B^ L — B <C Bp L . As the magnetic field approaches the critical value ( 1 1 
superconducting condensate becomes very small, 



the 



|(0)(B)|«tj, (12) 
and, consequently, the classical equation of motion (|7| can be linearized: 

{(® l -iD 2 )(D l +iD 2 )+e[B? L -B(x)}}cP=0. (13) 



The magnetic field inside the superconductor, B{x), appears explicitly in Eq. ( 13 i 
due to a rearrangement of the derivatives of the first term. 

In the vicinity of the critical magnetic field, B ~ B c , so that the second term in 
Eq. ( 13 i can be neglected and Eq. ( 13 1 is satisfied if 

£><£~0 with S=S 1 +/D 2 . (14) 



3.3.3 Lattice of Abrikosov vortices in background of magnetic field 

If the strength of the external magnetic field is smaller then the critical value ( [TT) 
then the superconductor may squeeze the magnetic field into the vortexlike struc- 
tures which are known as the Abrikosov vortices PT1 . The Abrikosov vortex is a 
topological stringlike solution to the classical equations of motion (|7]l and (18). A 
single Abrikosov vortex carries a quantized flux of the magnetic field, 
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d 2 x ± B(x ± ) = — , (15) 
e 

where the integral of the vortex magnetic field B is taken over the two-dimensional 
coordinates x 1 - = (x l ,x 2 ) of the plane which is transverse to the infinitely-long, strait 
and static vortex [notice, that the flux ( fl5] > is twice larger than a conventional value 
since we consider the condensed bosons <£> with the electric charge e and not 2e]. 

The Abrikosov vortex has a well-defined center where the condensate <P vanishes 
restoring the normal (nonsuperconducting) phase. At the vortex center the phase of 
the scalar field is singular. The behavior of the scalar field in the vicinity of the 
elementary vortex, situating at the origin and carrying the flux ( |15) , is as follows 

<P(x x ) oc \x- L \e ilf ' =xi+ix 2 , (16) 

where <p is the azimuthal angle in the transverse plane, and |jc x | is the distance from 



the vortex center. Equation dl6|l is valid provided m$|x x | -C 1 and m.j^\x | <C 1, 
where the mass parameters ma and m$ are given in Eq. d6j. 

If the external magnetic field is strong enough but still weaker than the critical 
value (IT) , then multiple elementary Abrikosov vortices may be created. Parallel 
Abrikosov vortices repel each other in a type-II superconductor, for which the mass 
of the scalar excitation is larger then the photon mass, w$ > m^. Due to the mutual 
repulsion, these vortices arrange themselves in a regular periodic structure known 
as the Abrikosov lattice [17|. The Abrikosov lattice corresponds to the so called 
"mixed state" of the conventional superconductor, in which both normal phase (in- 
side the vortex cores) and superconducting phase (outside the vortex cores) coexist. 

There are various arrangements of the Abrikosov vortices corresponding to dif- 
ferent types of Abrikosov lattices [17|. The simplest Abrikosov lattice is given by 



the following solution of Eq. ( 14 1 



<P = <P Q K(z/L B ) , K(z) = e -?<W 2 +; 2 ) £ e -^ 2 +2^ , (17) 



where <J>o is a dimensional complex parameter, z = x + ix , and 

L B = V2n£ B , l B = -^=, (18) 



is the inter-vortex distance L B expressed via the magnetic length l B . The area of an 
elementary lattice cell (i.e. of a cell which contains one Abrikosov vortex) is L B . In 



the solution ( 17 1 the vortices are located at the sites of the square lattice, 



^L=m+ 1 -, mez, i=i,2, (19) 

L B I 

at which the condensate <P{x\,X2) vanishes exactly. In the vicinity of these sites the 
scalar field (17 1 is described by Eq. ( fT6] ). 
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The ground state of the system in the mixed phase is characterized, by definition, 
by a minimal energy of the vortex lattice. One can show that a global minimum 
of the energy density of the GL model <|3j corresponds to a global minimum of a 
convenient dimensionless quantity which is called the Abrikosov ratio ifTTl : 



Ml 



(20) 



It turns out that the real ground state of the system corresponds an equilateral 
triangular lattice (which is sometimes also called "hexagonal" lattice) with the 
Abrikosov ratio /3a (Triangular) ps 1.1596, Fig. fright). For the square Abrikosov 
lattice the Abrikosov ratio ( 20 » is slightly higher, JSa (Square) pa 1.180, Fig.[3peft). 
Notice that despite very different visual appearances of these two lattices, the differ- 
ence in their energies (and in the corresponding Abrikosov ratios, /3a) is of the order 
of a few percent. A Ginzburg-Landau description of the type-II superconductors in 
the background of magnetic field can be found in the nice review |42). 



Lb 
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-3 ; 




1 

xi/L B 



-1 

Xi/L B 



Fig. 3 Minimal-energy vortex lattices for square (left) and equilateral triangular (right) lattices in 
the Ginzburg-Landau model. The dark dots correspond to the positions of the Abrikosov vortices. 
The equilateral triangular lattice corresponds to the global minimum of the energy. From Ref . 1 1 3 1 . 



3.3.4 London equations and complex electric conductivity 

The GL model ^ in the condensed phase describes a superconductor. Indeed, 
taking into account that in the ground state the condensate is a uniform time- 
independent quantity, one finds from from the definition of the electric current ([9]): 

d^-d v j£ L = -m 2 A F» v , (21) 



where niA is given in Eq. ([6|. Setting jj, =0 and v = 1,2,3, in Eq. (21 1 we recover 
the first London relation for a locally neutral [Jo(x) = 0] superconductor: 
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<?Jgl 
dt 



■ m A E, 



(22) 



where E or E' = —F°' is a weak external electric field. Equation ( 22 1 describes a lin- 
ear growth of the electric current in the constant electric field, thus indicating a van- 



ishing electric resistance of the system. According to Eq. ( 22 1, the superconductivity 
is homogeneous (coordinates-independent) and isotropic (direction-independent). 



The London equation ( 22 1 corresponds to a singular part of the complex conduc- 
tivity tensor Cu = Re Oki + 'Tm <%• The conductivity tensor is defined as follows: 



J k (x,t;co)= £otf(G))£/(x,*), 

k=l 



(23) 



where E(x, t) = Eoe'( x ' q fflr ' is the alternating external current in the long-wavelength 



limit, |q| — > 0. The London equation (22i indicates that (7u((o) — ol'" s ((Q 



OV., (ft)), where the singular part comes from the paired (superconducting) electrons, 



c» = 



nnr. 



8(co) 



2i 

KCO 



(24) 



while the regular part a M accounts for other contributions to the conductivity. 



3.3.5 Meissner effect 



The spatial components of Eq. pi) provide us with the second London relation: 

dxJ GL = -m 2 A B. (25) 
In the absence of a background electric field (E = 0), Eq. ^ implies Jgl = d x B, 

(26) 



so that Eq. (25 i can now be reformulated as follows: 

(-4+m|)B = 0. 



Equation (26 1 describes the Meissner effect: inside a superconductor the photon 
becomes a massive particle so that an external magnetic field B < B c is expelled. 
Physically, the Meissner effect appears because the external magnetic field induces 



circulating superconducting currents (25 i which, in turn, generate their own mag- 
netic field. As the generated field is directed in the opposite direction with respect to 
the external field, the magnetic field is eventually screened inside superconductor. 

If the external magnetic field is directed tangentially with respect the supercon- 
ductor's boundary then this field is always screened inside the bulk of the super- 
conductor. However, if the external magnetic field is directed along a normal of 
the boundary of a type-II superconductor, then the magnetic flux - which is a con- 
served quantity - may penetrate the superconductor and create a mixed phase of the 



Abrikosov vortices (Section 3.3.3 page 16 1. 
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3.4 Superconductivity of vacuum in strong magnetic field 



3.4.1 Electrodynamics of p mesons 

The conventional superconductivity is driven by the condensation of the Cooper 
pairs which are described by the local scalar field <P in the Ginzburg-Landau ap- 
proach. The superconductivity of vacuum in a sufficiently strong magnetic field is 
caused by emergence of quark-antiquark condensates which carry quantum num- 
bers of (charged) p mesons ifTUl . In order to describe the electrodynamics of the p 
mesons in the GL style, we use the following effective Lagrangian fl4l : 



1 p£V°* v +%?V 0) ' 1 - 



if = - \ V v - i(z^p v] )tz^.pv] +m 2 p y 

f—^pjl (27) 

where the complex vector field p^ = (p^ — r'p^)/-\/2 and the real- valued vector 

field p^ = p^ \ correspond, respectively, to the charged and neutral vector mesons 
made of the components of the triplet of the p field: 



Pf 



(2) (3; 
Pfi 



(28) 



The last term in Eq. ( 27 1 describes a nonminimal coupling of the p mesons to 
the electromagnetism via the field strength F^y = <3[^ A v ] of the photon field A^. 
The presence of the nonminimal coupling implies, in particular, the anomalously 
large v alue of the gyromagnetic ratio of the p meson, g = 2 (discussed already in 



Section 



3.1 



page 



13 I. Both the covariant derivative D u 



the strength tensor p^y = d^ p^ 1 — igsP^ Pv] involve the pnn coupling g s which 
has the known phenomenological value of g s ss 5.88, Ref. [ 14 1. 



,(0) 



(0) 



ieA^ and 



Lagrangian ( 27 1 is invariant under the electromagnetic gauge transformations, 



Jea(x) 



A^(x)+d^co(x), 



(29) 



which do not affect the neutral field, p^ ' (x) — > p^ ' (x) 



The p-meson Lagrangian (27i is an analogue of the Ginzburg-Landau La- 



grangian Q, while the p-meson field (28 i plays the role of the GL scalar field 
0. The electric current of the p mesons is given by the analogue of Eq. ([£]): 



Jp = /e[p v1 p VM - p v P VM +<5" \plpn~ pJiPv)] - -d v f, 



(0) 
Vfl ■ 



(30) 
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3.4.2 Instability of vacuum: potential energy in strong magnetic field 



The energy density £ of the p-meson ground state is given by the 7oo component, 

„2 



g = Too = 



1 



1 



(OK 2 



1 



2 - 0/ -r fry t ^ \H Ql ) ~*~ 4 (py° ) 2 + Y 



(pD 2 +(pD 



(0)\2 



t If 9 / + t \ e fOl e (0) 

+PoiPo,- + -pjjp,j+m z p {pipo+pjpt) ~ -foip'oi - 2 ^ I'aP,, ■ 



of the energy-momentum tensor of the p-meson electrodynamics (27 1: 



?uv — 2- 



HV ■ 



(31) 



(32) 



It is useful to consider a "homogeneous" approximation and ignore for a moment 
all derivatives and covariant derivatives in the energy density ( [3T) . This procedure 
corresponds, roughly speaking, to selection of a potential part of the energy density: 

1 2 4 

v(p M ,pf) = -B 2 + | £ [/(p / ip v -p^)] 2 + / e B(p 1 t p 2 -p 2 t pi) 



M.v=o 



,2 4 



(oh 



M=0 



E p^Pm • 



(33) 



M=0 



In the homogeneous approximation the ground state can be found via the minimiza- 
tion of the potential energy (33 i with respect to the fields p^ and pj±\ It turns out 
that in this approximation the vacuum expectation value of the neutral p-meson field 



is zero, p^' = 0. The quadratic part of the charged field is the following: 

2 

I 

1,6=1 



V (2) {p l i) = E Pl^abPb+m 2 p {plpo+plpi) 



m 2 p ieB 
-ieB m 2 p 



, (34) 



where the mass matrix for the Lorentz components p\ and p 2 is non-diagonal. 

The eigenvalues fi± and the corresponding eigenvectors p± of the mass matrix 
(34 1 are, respectively, as follows: 



Hi = m 2 p ±eB, p± = -^= (pi T ipi) 



1 



(35) 



It is clearly seen that one of mass states, either /i_ or ji + depending on the sign of 
eB, is getting smaller as the magnetic field increases. Taking for definiteness eB > 0, 
we chose the ground state of the system in the following form: 



pi=p, p 2 = -/p, Po = 0, p 3 =0. (36) 
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The longitudinal components po and P3 are always zero because for any value of 



the magnetic field the corresponding terms in Eq. ( 34 1 are positive and diagonal. 
The total potential energy of the p meson system in the homogeneous ground 



state can be calculated with the help of Eqs. ( 33 i and ( 36 1 



v(p) = \b 2 



where 



-2e{B c -B)\ P \ 2 +2g 2 s \ P \ A 



(37) 



B c 



(38) 



is the critical magnetic field ([TJ. 

Thus, we get the familiar Mexican-hat potential ( 37 1 for the p-meson conden- 
sate p. In particular, the very same form of the potential appears in the GL model 
of superconductivity ( 10 1, with one very important exception: in the conventional 



superconductivity the condensation of the Cooper pairs emerges at low tempera- 
tures, T < T c , while the p mesons start to condense in the vacuum in the presence 
of sufficiently strong magnetic fields, B > B c : 



I(P)K 




e(B-B c ) 



B<B C . 



(39) 



Here the subscript V indicates that we consider the potential part V(<P) only so that 
we ignore all kinetic terms. If B > B c then the condensate ([39} breaks spontaneously 



the electromagnetic symmetry group ( 29 1, similarly to the spontaneous symmetry 



breaking in the superconducting phase of the GL model. 



3.4.3 Negative condensation energy due to p-meson condensate 



The homogeneous nature of the p-meson condensate (39 1 is an artifact of the po- 
tential approximation which was useful to find the very presence of the tachyonic 
instability of the noncondensed state at B > B c . In order to determine a detailed 
structure of the ground state, we should work beyond the potential approximation. 
To this end, we notice that a wavefunction of the lowest energy state of a free parti- 
cle in a uniform static magnetic field is independent on the coordinate z = x 3 along 
the magnetic field axis. Secondly, the dependence on the time coordinate t = x° 
should appear only in a form of a trivial phase factor. Thus, we concentrate on x 1 - 
and x 2 -dependent solutions to the classical equation of motions for the p mesons, 
similarly to the case of the ordinary Abrikosov lattice solutions. 

Technically, it is convenient to choose the complex coordinate z — x 1 + ix 2 and 
define the complex variables 6 = G\ + 162 an d their conjugates 6 = G\ — Wi 
for vector quantities 0^. Then, the classical equations of motion of the p-meson 



model ( 27 1 can be written in the following "complexified" form: 
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gs dB + iemlp {0) =0, (40) 
{-dd+m 2 + 2g 2 s \p\ 2 )p {0) -2ig s d\p\ 2 -0, (41) 
[-DD + 2(g s C-eB + 2g 2 \p\ 2 +m 2 p )]p=0, (42) 

where the z-components of the magnetic field and its analogue for the neutral p 
mesons are, respectively, as follows 

B{z) = d,A 2 - d 2 A, = Im(dA) , C(z) = <?ipf ' - d 2 pf ] = lm(dp^) . (43) 
We have also introduced two covariant derivatives: 

D=D 1 +iD 2 =T) + ig s p i0) , T) = d-ieA, A = ^±z, (44) 

2i 

where the external uniform magnetic field Z? ext should be distinguished from the 
full magnetic field in the superconducting state, B(x\ ,x 2 ) = B(x\ + ix 2 ) as the latter 
induces a backreaction from the superconducting ground state. 

The transverse components of the electric current p0|) are as follows: 



7 X =/i +U 2 =2/e(p t Dp + <9|p| 2 ) +/— dC. (45) 

gs 

In the vicinity of the phase transition, B e xt > B c with \B ext — B c \ B c , the equa- 
tions of motion ( [40] ), ( pT) and ( |42| can be linearized. It turns out that the equation for 
the p -meson condensate in the overcritical magnetic field (B ext > B c ) in the vacuum 
is identical to the equation for the Cooper pair condensate ( 14 1 in the subcritical 
magnetic field (B ext < B c ) in the GL model of conventional superconductivity 1 10 1: 

Dp = (d- e -B ext z)p=0. (46) 



Among infinite number of solutions to Eq. (46 1, the ground state solution corre- 
sponds to the global energy minimum. In the chosen approximation, the mean value 
of the energy density (31 1 can be expressed via the p-meson condensate fl3l : 

{S) ^ (Too) = l -Bl xt +2e(B c -B eM )(\p\ 2 )+2e 2 (\p\ 2 ) 2 

+2(g 2 -e 2 )(\p\ 2 -^- 2 \p\ 2 ), (47) 

where d 2 = d 2 + d 2 is the two-dimensional Laplace operator, 

1 , 1 _ . . , 



(x ± ) = — XoKKI) (48) 



-d\+ml v ~ ' 2% 



is a two-dimensional Euclidean propagator of a scalar massive particle with the mass 
of the neutral p '°' meson, 
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ra 



: m p(0) 



„2v _I 
2 



(49) 



and Kq is a modified Bessel function. Contrary to the potential part of the energy 



density ( 37 1, the full expression (47 1 depends nonlocally on the condensate p. 



A general solution of Eq. ( 46 1 is the following generalization to the square 
Abrikosov lattice ( fT7| i: 



P(z)= E C„exp|--(|; 



2 2 

7TV n 



z = xi + «2 , (50) 



where v is an arbitrary real parameter, Lb is the magnetic length ( 18 1 and C„ are 
arbitrary complex coefficients. 

The ground state solution is given by an equilateral triangular lattice, Fig. fright). 
The corresponding coefficients C„ obey the two-fold symmetry C„+2 = C„ with 
C\ = /Co, while the independent parameters v and Co cannot be calculated analyt- 
icity so that they should be found by a numerical minimization of the energy den- 



sity ( 47 1. Equivalently, one can also minimize an analogue the Abrikosov ratio ( 20 1 



|P| Z 
(IPl 2 ) 



The left and right panels of Fig.|4]show, respectively, the condensation energy 



(51) 



8S = (<?) 



(52) 



and the superconducting condensate \p \ = sj (|p| 2 ) in the ground state. The rise of 
the condensate at B > B c makes the ground state energy smaller compared to the 
normal, noncondensed state. 




B/Bc B/B c 



Fig. 4 At B > B c the superconducting state - corresponding to the equilateral triangular lat- 
tice Fig. [5Jright) - is more energetically favorable compared to the trivial vacuum state: at the 
critical magnetic field B = B c , the condensation energy {52} , the left panel, becomes negative due 
to emergence of the superconducting condensate (the right panel). From Ref. 1131 . 
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3.4.4 Periodic pattern of ground state: superfluid and superconductor vortices 



Solutions of Eq. ( 50 » are inhomogeneous functions in the transversal {x\,X2) plane. 
The inhomogeneities in the charged p condensate induce an unexpected condensa- 
tion of the neutral p mesons: 



d\p\ 2 ^d J ^V^o(m |^-y x |) |p(/)| 2 . (53) 



The longitudinal components of the neutral condensate are zero, pl°^ = pi ^ = 0. 

The external magnetic field B e xt leads to a backreaction from the charged con- 
densate d50), which creates a transverse electric current d45b, 



mo\x ± -y ± 



2iem} x d 

and, consequently, affects the magnetic field inside the superconductor: 



,,,, / Hemkd , ,,\ . , . iernk _ f , , 



|P(/)| 2 (54) 



B(x ± )=B ext + - 



p{x^\ 2 -{\p\ l ) 



(55) 




Fig. 5 The charged superconducting (left) and neutral superfluid (right) condensates in the 
transversal (x\,X2) plane at B = 1.01B C . The 3D plots illustrate the absolute values of the con- 
densates (in MeV) while the corresponding projections on the {x\,x%) planes of these figures are 
the density plots of the phases of the corresponding condensates. The white lines of the projections 
are gauge-dependent singularities (the Dirac sheets) which are (left) attached to the superconductor 
vortices and (right) stretched between the superfluid vortices and antivortices. 

In Fig. [5] the charged and neutral condensates are plotted as functions of the 
transverse coordinates x\ and X2 for the magnetic field B = 1.01B C . The periodic 
equilateral-triangle structure of the absolute value of the charged p-meson conden- 
sate is identical - apart from the magnitude of the physical scales - to the one of the 
GL model, Fig. fright). The absolute value of the neutral condensate also exhibits 
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a lattice pattern which has, however, a bit more involved appearance: hexagonally- 
shaped structures are arranged into the equilateral triangular lattice, Fig. [6] 



Fig. 6 A visualization of the 
nested structure of the electri- 
cally charged, superconduct- 
ing condensate \50\ and the 
electrically neutral, superfluid 
condensate l |53[ l, plotted in 
magenta and green, respec- 
tively. Shown are the regions 
where these condensates take 
maximal values. The trans- 
verse plane corresponds to 
2fm x 2fm region at the 
magnetic field B = 1.01B C . 





The charged and the neutral p-meson condensates coexist together. Since the 
magnetic field cannot directly induce the neutral condensate, the mechanism of its 
appearance is as follows: the background magnetic field induces the charged con- 
densate p, Eq. (50 1, while the charged condensate gives rise to the neutral one, p(°\ 



Eq. (53 I. As a result, the neutral condensate is an order of magnitude smaller than 
the charged condensate, Fig. [5] These condensates form a nested structure, Fig. [6] 

The projection of Fig. |5](left) shows the density plot of the phase of the charged 
condensate, argp. The phase - which is not a periodic function of the transverse 
coordinates X\ and x 2 - exhibits discontinuities across which the phase is changed 
by 2k. These discontinuities correspond the Dirac sheets, shown as the white lines in 
the same projection of Fig.|5](left). The Dirac sheets are attached to the new class of 
vortices, "the superconductor p vortices". The positions of these vortices correspond 
to the endpoints of the Dirac sheets. According to the 3D plot of the same figure, 
the absolute value of the p-meson condensate is vanishing at the centers of the 
superconductor vortices. Locally, these superconductor vortices have the structure 
which is similar (up to a gauge-dependent phase) to the Abrikosov vortices in the 

: p (x 1 - ) = pi (x x ) = ip 2 (x 1 - ) oc \x x \ ( 



16) 



\x\ -(-un- 



conventional superconductors 

Contrary to the phase of the charged condensate pO) , the phase of the neutral 
condensate ( |53] > is a periodic function of the x\ and x 2 coordinates. The neutral 
phase exhibits the 27T-discontinuities as well. These discontinuities are visualized 
as white lines in the projection on the bottom-left panel of Fig. [5] The end-points 
of these discontinuities mark positions of a new type of vortices called "superfluid 
p vortices" connected by the corresponding 27T-discontinuity to the superfluid an- 
tivortices. Locally, the superfluid vortices have, up to a phase, the familiar structure: 
p (°) (x L ) = p| 0) (x x ) + ipf ) (x 1 - ) °c Ix 1 - \e i <P=xi+ ix 2 . 

The vacuum ground state has a rich "kaleidoscopic" structure in terms of the 
vortex content: the equilateral triangular lattice of the superconductor vortices is 
superimposed on the hexagonal lattice of the superfluid vortices and antivortices. 
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Fig. 7 (left) The kaleidoscopic vortex structure of the charged and neutral condensates induced 
by the magnetic field B = 1.01S C : the superconductor vortices (the large red circles) are always 
superimposed on the superfluid vortices (the small blue disks marked by the plus signs) forming an 
equilateral triangular lattice in the (x\,x%) plane. The isolated superfluid vortices and antivortices 
(the small yellow disks with the minus signs) are arranged in the hexagonal lattice. The shades of 
green illustrate the absolute value of the neutral p-meson condensate {53} (from Ref. 1 13 1). (right) 
The density and the vector flow of the superconducting currents in the (x\,X2) plane at B = 1.01 B c . 



An elementary lattice cell of the kaleidoscopic lattice contains one superconductor 
vortex in the electrically charged p condensate as well as three superfluid vortices 
and three superfluid antivortices in the neutral p(°) condensate, Fig.^left). 



3.4.5 Superconductivity and superfluidity in the ground state 

Electric transport properties of a material (such as, for example, the electrical con- 
ductivity) are usually determined in a linear response approximation in which one 
studies an electric current generated inside the material by a weak (test) external 
electric field background. The electric field should be weak enough in order to pre- 
serve, in a leading order, the ground state of the studied material. 

In our ground state the transverse (with respect to the strong magnetic field) 



electric currents ( 54 1 of charged condensates are confined to elementary cells of the 
periodic ground state, Fig. [TJright). The size of the elementary cell is of the order 
of the size of the wavefunction of lowest Landau levefl In order for a net electric 
current to be induced in the transverse directions, the quarks need to be excited to 
the next Landau level which is, however, separated from the lowest Landau level by 
a large energy gap of the order of 8E ~ v/|eB|. It is the energy gap which makes the 
vacuum state to behave as an insulator in the transverse directions because a weak 
(|E| <C |B|) transverse (E _L B) electric field E cannot create large enough excitation 
of overcome the gap. The presence of the gap is the very reason why the Meissner 



7 In physical units the size of the cell is approximately 0.5 frn for the near-critical field, B ~ B c 
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effect is absent in the superconducting ground state [ 10 1 so that the emerging super- 



conductivity does not screen the external magnetic field (Section 2.2.3 pagep 





Fig. 8 The strength of the vacuum (left) superconductivity K, Eq. |57| , and (right) superfluidity 
k' ', Eq. |59| l, as the function of the transverse coordinates x\ andx2 at magnetic field B = 1.01 B c . 
In addition, the left plot illustrates the superconductor vortices (the large red tubes) and the super- 
fluid vortices and antivortices (the smaller blue and yellow tubes, simultaneously) in accordance 
with notations of Fig^left). In the right plot the semitransparent plane highlights the line k' ' = 
where the superfluid strength changes its sign. 



Contrary to the transverse electric currents, the longitudinal currents are not 
restricted by the external magnetic field. Let us apply a weak electric field E = 
(0,0, £V) along the axis of the strong magnetic field B = (0,0, B). According to the 



equations of motion of the p-meson model ( 27 1, the induced electric currents ( 30 1 
satisfy the following equations 



djj,(x) dJo(x) 



-k(x ± )E z 



dJ k {x) _ dJy (x) 
dxV dx k 



•0, 



(56) 



where \i = 0, . . . , 3 and k = 1,2. 



The set of equations (56 1 is nothing but an anisotropic "vacuum" analogue of the 



London equation pi 
current 



of superconductivity. Equations (56i show that the electric 
field 



induced by a weak electric "test" field - flows without resistance along the 
magnetic field axis while in the transverse directions the superconductivity is absent. 
The strength of the vacuum superconductivity is characterized by the quantity K, 
which is a nonlocal function of the superconducting p-meson condensate: 



(* X ) = (- 



K(X 



4e 2 ml 
-d 2 +m 



(x 1 



K (. 



TOO U 



-y 



p(y X )| 2 -(57) 



According to Fig. [FJleft), the strength of the superconductivity (57 1 is a weakly 
dependent function of the transverse coordinates x\ and xj_. In a response to a weak 
electric current, the superconducting currents are generated outside of the supercon- 
ductor vortex cores and the maxima of the induced electric currents are located at 
the centers of the superfluid vortices. Contrary to the ordinary superconductivity, 
the vacuum superconductivity is not completely suppressed inside the vortices due 
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to the nonlocal nature of the relation between the transport coefficient K and the 



superconducting condensate ( 57 i 



Unexpectedly, the condensate of the neutral p mesons is also sensitive to the 
presence of the external electric current. It turns out that the electrically neutral 
current of the p' ' mesons, defined via the relation jjy = — |-<3 v /v/i, satisfies a 
London-like equation as well B31I : 



kW(x u x 2 )E z (x), (58) 



dx° dx 3 

K®(x lt xi) = ( A £ d } 2 \p\ 2 )(x u x 2 )^^dlK(x^, (59) 

where the superfluid coefficient is visualized in Fig. fright). 

Equations ( |56| and ( |58] l indicate that, respectively, the charged and neutral cur- 
rents should flow frictionlessly (i.e., accelerating ballistically) along the magnetic 
field axis if an weak external electric field is applied along the same direction. Then, 
if even at some moment of time the electric field is set back to zero, both the su- 
perconducting current and the superfluid flow would continue to flow permanently 
because of the absence of the dissipation forces for the corresponding condensates. 

Notice that the electric-field-induced superfluid flow is a locally nonvanishing 
quantity, while the total superfluid flow of each elementary lattice cell is zero, 

<?(/f)xM d(4 0) )x(*)_ n 



fjv „ (jv ;, =0, {ff) x = Jd^ff{x), (60) 

because the cell-averaged superfluid coefficient \59) is zero, {k ( -°\x j -)} ± =0, too. 
A comparison of the left and right plots of Fig 7|8] reveals that the external elec- 
tric field generates the positive superfluid flow at the positions of the superconduc- 
tor vortices which are always accompanied by a superfluid vortices according to 
Fig. [Tpeft). The negative superfluid flow is generated at the positions of other, un- 
accompanied superfluid vortices and antivortices. 

Finally, we would like to stress that the global quantum numbers of the new 
exotic superconducting (and, simultaneously, superfluid) phase correspond to the 
quantum numbers of vacuum. For example, all chemical potentials in the super- 
conducting phase are vanishing. The vacuum is an electrically neutral object: the 
presence of the positively charged condensate p implies an automatic appearance of 
a negatively charged condensate p * of the equal magnitude, p = |p * | . As a result, in 
strong magnetic field the energy of the vacuum is lowered due to the emergence of 
the charged condensates, while the net electric charge of the vacuum stays always 
zero iTTOl fTTTl . Despite of the net electric neutrality, the vacuum exhibits the super- 
conductivity since a weak external electric field pushes the positively and negatively 
charged condensates in opposite directions along the magnetic field axis, thus cre- 
ating a net electric current of a double magnitude. 
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3.4.6 Abelian gauge symmetry breaking and gauge-Lorentz locking 

What is the symmetry breaking pattern in the new superconducting phase of the vac- 
uum? The vacuum superconductivity appears due to the emergence of the magnetic- 



field-induced p-meson condensate (36 1. In the presence of the background mag- 



netic field B, the group of global rotations SO(3) tot of the coordinate space is ex- 
plicitly broken to its <9(2) lot subgroup which is generated by rotations around the 



magnetic field axis. The scalar combination p of the vector condensates (36 1 trans- 
forms under the residual rotational group 0(2) lot as follows: 

0(2) rot : p(x)->e?*p(x), (61) 

where <p is the azimuthal angle of the rotation in the transverse plane about the xj 



axis. The p meson field transforms also under the electromagnetic group (29 1: 

t/(l)e.m.: pW^^lpW, (62) 

Thus, if the condensate p were a homogeneous (i.e., coordinate independent) quan- 
tity then the ground state would "lock" the electromagnetic gauge symmetry with 
the residual rotational symmetrjj^] t/(l) e . m . x 0(2) ml — > I/(l)i cked> since a rotation 
of the coordinate space at the angle <p about the axis xj and a simultaneous gauge 
transformation with a constant gauge-angle co — —(p/e leave the homogeneous con- 
densate p intact. The inhomogeneities in the p condensate break the locked sub- 
group further from the global 1/(1) group down to a discrete subgroup of the lattice 
rotations G"}^^^ of the kaleidoscopic lattice state, Fig. gleft), Ref. IT51 : 

f/(l) e . m .xO(2) rot ^G^ cked . (63) 



3.5 Superconductivity of vacuum in Nambu-Jona-Lasinio model 

Basic properties of ordinary superconductivity can equally be revealed either in 
the Ginzburg-Landau (GL) phenomenological approach which describes the scalar 
field of the superconducting carriers or in the Bardeen-Cooper-Schrieffer (BCS) 
model which accounts the electron pairing into the superconducting carriers (Sec- 
tion : 



3.2 page 14 1. Both approaches are mathematically equal if the temperature is 
sufficiently close to the superconducting phase transition ll33l . 

So far we have discussed the basic features of the vacuum superconductiv- 
ity in the effective p meson electrodynamics lfl4l . which serves as a "vacuum" 
analogue of the GL approach to the ordinary superconductivity (34). Below we 
briefly consider, following Ref. ifTTI . the p-meson condensation in the BCS-like 
approach 11371 . which is based on the Nambu-Jona-Lasinio model l36l[T5l . 



8 A philosophically similar phenomenon, a color-flavor locking, is realized in a different context 
of the color superconductivity in a dense quark matter 1 44 1 . 
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3.5.1 Effective action in strong magnetic field 

We consider an extended two-flavor (Nf — 2) three-color (N c — 3) Nambu-Jona- 
Lasinio (NJL) model [15|: 



4^3 

2 



(64) 



where the light quarks are represented by the doublet y = (u,d) T , and G s and G v 
are the corresponding bare couplings of scalar and vector quarks' interactions. The 
masses m u and m^, and electric charges (q u = +2e/3 and q ( / = —e/3) of the quarks 
are combined into the bare mass matrix M° = diag(m°,mji) and the charge matrix 
Q = d.ia.g(q u ,q c i), respectively. The 2x2 matrices in the flavor space are denoted by 
hats over the corresponding symbols and T ! , i = 1,2,3, are the Pauli matrices. 
The partition function of the NJL model can be represented as an integral, 

2T = f DyDy e ^ d " xje = [ DaDnDVDAe^^ , (65) 



over bosonic fields [ 15 1. The bosonic fields are given by one scalar field a ~ yy, 
the triplet of three pseudoscalar fields % ~ yfxy [made of the electrically neutral, 
7t° = 7T 3 , and electrically charged, ft* = (n l =p i7t 2 )/V2, pions], four vector fields, 

[composed of the flavor-singlet coordinate-vector 0)-meson field co^, and of the 
electrically neutral, = p 3 , and charged, p* = (p^ =p ip^)/y2, components of 
the p -meson triplet], and four pseudo vector (axial) fields, 

where the fields and represent the singlet axial f\ meson and the ai 

triplet of the axial mesons, respectively. 

The effective bosonic action in Eq. (65 i is as follows 

(68) 



S[a, k, V,A] = S„ - [d 4 x l-Lr (a 2 + % 2 ) * (V^ +A k ^) 



2Gf V ' 2Gf 



S V [<J, tz,V,A] = -iN c TThn(m) , (69) 
iD = + Q $ - M° + V p + f4 - (<r + if tct) . (70) 
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where we have used simplified notations for the expectation values of the fields, 
(a) = C etc. In the absence of a magnetic field background the expectation value 
of a plays a role of the constituent quark mass, m q = a ~ 300 MeV while the 
expectation values of the fields n, V, and A are zero fl5l . 



The effective action ( 68 1 in the strong magnetic field background was calculated 
in Ref. [11] in the lowest Landau level (LLL) approach using a mean-field technique 
inspired by calculations of the magnetic catalysis phenomenon [1J. In the regime of 
the LLL dominance the propagator of a /'s quark 



S^ L (x,y)=Pjr^y)s\(J-y^), 
factorizes into the B-transverse projector onto the LLL states 

P ± (x ± V X ) - kMAlf^ab^-i^fB^-y^) 2 

and B-longitudinal fermion propagator in the 1 + 1 dimensions, 

D ii ^-iffy 2 



4(*||) = in(*l|) 



1\ Pf ■ 

yll h\ —m J 



f 



(71) 



(72) 



(73) 



which depend, respectively, on the B-transverse, x x = (x 1 ,x 2 ), and Z?-longitudinal, 
x" = (x°,x 3 ), coordinates [lj. Here qf is the electric charge of the / th quark (/ = u 1 d) 
and eB > is ta ken for definiteness. 



In Equation (73 i the matrix Pj is the spin projector operator onto the fermion 



states with the spin polarized along (for u quarks) or opposite (for d quarks) to the 
magnetic field (we use / = ±1 for, respectively, / = u,d). The operator pj projects 
the original four 3+1 fermionic states onto two (l+l)-dimensional fermionic states, 
so that in the LLL approximation the quarks can move only along the axis of the 
magnetic field (the latter fact is a natural sequence of the LLL dominance IT]). 



The operator (72 1 satisfies the projector relation, 



P f oP f= P f> 



AoB 



d\ ,\i....v )«(v ....). 



(74) 



where "o" is the convolution operator in the Z?-transverse space. 

In the one-loop order the effective action (|69| contains a scalar and vector parts: 



(75) 



S = S s (o,7l)+S v (A,V), 



respectively. In terms of the nontrivial condensate^] the potential term in the scalar 
part of the action has the following (renormalized) form: 



Ss = - 



d 4 x 



1 

2Gs ( 



\eB\Nc 

87T 2 



In ^ - 1 ) a 



(76) 



Here we omit all terms with vanishing condensates as well as all kinetic terms. 



Electromagnetic superconductivity of vacuum induced by strong magnetic field 



33 



which reflects one of the most important features of the magnetic catalysis [ 1 J: an 
enhancement of quarks' masses by the magnetic field background, 



m q {B) = a mm (B) ^^exp{-27Z 2 /{G s N c \eB\)} , 



given by the minimum a m j n of potential (76i. The mass scale, fi °c 



(77) 



is to be 



fixed beyond the LLL approximation because it is determined, in particular, by the 
(1 + 1) dimensional motion of the quarks along the magnetic field Qj. As noticed 
in Ref. [ 1 1, the renormalization of the scalar NJL coupling Gs in the MS scheme, 



1 



1 N c \eB\ 



G s 



27tGgn 
N c \eB\ 



(78) 



is very similar to the renormalization of the coupling constant Ggn in the 1+1 
dimensional Gross-Neveu model |38|. The divergencies of the 1 + 1 dimensional 
fermions are treated in the dimensional regularization in d = 2 — 2e dimensions, 
l/e = 1/e — Ye + log47z: and y E s» 0.57722 is Euler's constant. 

A potentially nontrivial part of the (non-renormalized yet) effective vector action, 



-,(0) _ iNc 



Tr 



IX>0 l£>0 



1 a' 

= — In — 
e p z 



(<j>*oP e o(j))- 



AN c \eB\ 
9n 2 

,2 



(79) 



l-ln^-2)(rof'.o<) 



involves only the B-transverse combinations of the vector and axial mesons, = 



(Pj + ipj )/2 and | = (a \ + iat)/2. In Eq. (79 1 the B-transverse projector for the 



unit charged particle P e is given by Eq. (72 1 with the replacement qj — > e: 

9k 



eB 



(80) 



The potential ( 79 1 has an unstable tachyonic mode which is determined by an 
inhomogeneous eigenstate of the charge-1 projection operator ([80|: 



( j P e o0)(x ± ) = 0(x ± ). 



(81) 



The solution to this equation is a general periodic Abrikosov-like configuration ifTTl 



which is given, up to a proportionality coefficient, by Eq. ( 50 »: <j>(x ) « p(x ). One 



can also show [11] that there are no unstable modes for the axial mesons and, in 
accordance with Eq. ( [36| , no unstable modes exist for B-longitudinal components 
of the p mesons. Thus, we set the corresponding expectation values to zero. 

For the sake of simplicity, we set in Eq. ( 50 1 all coefficients C„ equal, C„ = <po, and 



then we evaluate certain basic quantities for the simplest square lattice ( 17 1. As we 



have mentioned, despite different visual appearances of the square and equilateral 
triangular lattices, Fig. [3] the corresponding bulk quantities (as, for example, the 
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energy density) evaluated at these condensate solutions are almost the same as the 
difference between them lies within (sometimes, a fraction of) percents. 

The leading quadratic and quartic terms in the effective potential for the square 



lattice solution ( 17 1 of the eigenvalue equation dSTb are given by 



V=V2 



1 

G~b 



2N c \eB\ 
9n 2 



\<h\ 



r |4 



1 



1 



9G S 



(82) 



where Co ~ 1.2 is a numerical (geometrical) factor. If the magnetic field exceeds 
certain critical strength^] eB™ L = 9n 2 /(2N c G B ) ~ 1 GeV 2 , then the potential (82 1 
becomes unstable towards a spontaneous creation of the B-transverse p-meson con- 
densates with the tachyonic mode pj~ (x^) — —ip^ ( x± ) — 0( x± ) [and , respectively, 



Pi(x ) — ip^ (x ) = <j)*(x )], where <j)(x ) is a solution of Eq. (81 



3.5.2 Electromagnetically superconducting ground state in the NJL model 



In the magnetic field background, the effective p-meson potential in the NJL 
model d82l has the same Ginzburg-Landau form as the potential ( 37 1 for the p- 



meson field in the p-meson electrodynamics (27 I. If the magnetic field exceeds the 
critical value, B ^ B^ 1L , then the charged p-meson condensate emerge. In terms of 
the quark fields the new vector quark-antiquark condensates are as follows: 



(ujid) — ~i{uj2d) 



MB) 



K 



tx\ +ix2\ 



(83) 



where the function K(z) is given in Eq. ( 17 I. The magnitude and global phase 8q of 



the condensate ([83[l are determined by the following formula: 



MB) = e l0 °C^m q {B) 1 - 



B 



NJL\ 1/2 



B 



for B > B 



NJL 



(84) 



Here Ca_~ 0.5 1 is a numerical prefactor and the /^-dependent quark mass m q is given 
in Eq. ( 77 1. At B < B^ 1L the p-meson condensate ( 84 1 is zero. The superconducting 



phase transition at B — B c is of the second order with the critical exponent 1 /2, 



similarly to the phase transition in the p-meson electrodynamics (27 I. 

The quark condensates ( |83"] l have the quantum numbers of the p mesons. They 
form an inhomogeneous ground state identical to the one found in the p-meson 



electrodynamics (50 1. It is very interesting to notice that the ground state in the NJL 
model (64 1, determined by the integral equation ( 81 1, and the ground state in the p- 



meson electrodynamics (27 1, determined by the differential equation of motion (46 1, 
have exactly the same functional form d50]l. 



We have estimated the critical field only approximately since the phenomenological values of 
the NJL parameters Gsy are not known precisely |45 1. Moreover, subtleties of the renormalization 
of the effective dimensionally reduced (1 + 1) -dimensional theory embedded in 3+1 dimensions 
provide us with an additional uncertainty. 
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The vacuum state (|83]l of the NJL model is superconducting. The validity of the 



anisotropic London equation (56 1 for the quarks' electric current, 

/M W = Lf^MW^Vf) = -tttfQSfrx)] , (85) 
can be shown in a linear-response approach using retarded Green functions ifTTl : 



d(J 3 ){J) d(J )(x\\ 



<?x 3 



2C q 
(2n)- 



e 3 (B-B™ L )E 3 for B > B, 



NJL 



(86) 



where C„ ~ 1 is a numerical prefactor and at B < B c the right hand side of Eq. ( 86 1 



is zero. For the sake of simplicity, in Eq. ( 86 1 we have averaged the electric charge 
density J° and the electric current J z = J° over the transverse plane ( 60 1. 



Apart from the prefactors, the London equations in the NJL model ( 86 1 and in 



the p-meson electrodynamics (56 1 are identical. In a linear-response approximation 
these laws can be generalized to a completely Lorentz-covariant form ifTTIl . 



dfiJ v - dvJ^ = y- (F,F)F^ V , 



(87) 



via the Lorentz invariants (F,F) = 4(B,E) and (F,F) = 2(B 2 -E 2 ). Here y is 
a function of the scalar invariant (F,F) and F^ v = e^yapF"^ /2. The Lorentz- 
covariant forms of the local London laws for the superconductor ( 56 1 and super- 
fluid (|58]l components can be rewritten in a similar way. 



4 Conclusion 

We have shown that in sufficiently strong magnetic field the empty space becomes 
an electromagnetic superconductor. The new state of the vacuum has many unusual 
features lfT0lfTT1[T3ll4ll: 

• The magnetic field induces the superconductivity instead of destroying it. 

• The Meissner effect is absent. 

• The superconductivity has a strong anisotropy: the electric currents may flow 
without resistance only along the axis of the magnetic field. 

• The superconductivity appears in the empty space as a result of the restructuring 
of the quantum fluctuations due to the presence magnetic field. The overcritical 
magnetic field (B > B c ~ 10 16 T) induces the quark-antiquark condensates which 
have the quantum numbers of the electrically charged p-mesons. 

• The electromagnetic superconductivity is always accompanied by the superfluid 
component caused by emergence of a neutral p-meson condensate. 

• The tandem superconductor-superfluid ground state is inhomogeneous, it resem- 
bles an Abrikosov lattice in a mixed state of an ordinary type-II superconductor. 

• The charged and neutral vector quark-antiquark condensates have stringlike topo- 
logical singularities: superconductor and superfluid p vortices, respectively. 
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• The ground state is characterized by a "kaleidoscopic" lattice structure made of 
the equilateral triangular lattice of the superconductor vortices which is superim- 
posed on the hexagonal lattice of the superfiuid vortices. 

The vacuum superconductivity may be considered as a "magnetic" analogue of 
the Schwinger effect. Indeed, the Schwinger effect (the vacuum superconductivity) 
is the electron-positron pair production (the emergence of the quark- antiquark con- 
densates) due to strong electric (magnetic) field background in the vacuum. Con- 
trary to the Schwinger effect, the vacuum superconductivity is a state, not a process. 

The sufficiently strong magnetic fields, of the strength from two to three times 
higher than the required critical value ([T| may emerge in the ultraperipheral heavy- 
ion collisions at the Large Hadron Collider (LHC) at CERN [5 1. Thus, signatures of 
the magnetic-field-induced superconductivity have a chance to be found in labora- 
tory conditions. 
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